Moments of nucleon generalized parton distributions from lattice QCD 
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We present results on the lower moments of the nucleon generalized parton distributions within 
lattice QCD using two dynamical flavors of degenerate twisted mass fermions. Our simulations 
are performed on lattices with three different values of the lattice spacings, namely a = 0.089 fm, 
a = 0.070 fm and a = 0.056 fm, allowing the investigation of cut-off effects. The volume dependence 
is examined using simulations on two lattices of spatial length L = 2.1 fm and L = 2.8 fm. The 
simulations span pion masses in the range of 260-470 MeV. Our results are renormalized non- 
perturbatively and the values are given in the MS scheme at a scale /i = 2 GeV. They are chirally 
extrapolated to the physical point in order to compare with experiment. The consequences of these 
results on the spin carried by the quarks in the nucleon are investigated. 



PACS numbers; ll.15.Ha, 12.38.Gc, 12.38.Aw, 12.38.-t, 14.70.Dj 



I. INTRODUCTION 



Lattice QCD calculations of observables, related to the structure of baryons, are now carried out using simulations 
of the theory with parameters that are close enough to their physical values that a connection of lattice results to 
experiment is facilitated. This is due to the fact that systematic uncertainties caused by a finite volume, a finite 
lattice cut-off and unphysically high pion masses are becoming better controlled. Nowadays, a number of major 
collaborations are producing results on nucleon form factors and the first moments of structure functions close to the 
physical regime both in terms of pion mass and with respect to the continuum limit [H-Q . 

The Generalized Parton Distributions (GPDs) encode important information related to baryon structure f^^. 
They occur in several physical processes such as Deeply Virtual Compton Scattering and Deeply Virtual Meson 
Production. Their forward limit coincides with the usual parton distributions and their first moments are related to 
the nucleon elastic form factors. Moreover a combination of their second moments, known as Ji's sum rule [lo| . allows 
to determine the contribution of a specific parton to the nucleon aiigular momentum. In the context of the "proton 
spin puzzle" this has triggered an intense experimental activity [Ill4l5| . 

The GPDs can be accessed in high energy processes where QCD factorization applies, in which case the amplitude 
is the convolution of a hard perturbative kernel with the GPDs, as illustrated in Fig. [TJ Generically the GPDs are 
defined as matrix elements of bi-local operators separated by a light-like interval. Due to the Wick rotation such 
matrix elements cannot be computed directly on the Euclidean lattice. Instead one considers their Mellin moments, 
which in principle, carry the same information. 

If \p') and \p) are one-particle states, the twist-2 GPDs, which are studied in this paper, are defined by the matrix 
element [T6j : 

A/2 

1 f rl\ — ^9 f dan-A{na) 

Fv{x,i,q^)^-j—e"''^{p'\i^{-\n/2)TVe-^'^ V'(An/2)b) , (1) 
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where q = p' — p, C = '9/2, x is the momentum fraction, and n is a hght-hke vector colhnear to P = {p + p')/2 and 
such that P ■ n — 1. The gauge hnk 7-'exp(. . . ) is necessary for gauge invariance. In model calculations it is often set 
to one, which amounts to working with QCD in the light-like gauge A ■ n = 0, but on the lattice such a gauge fixing 
is not necessary. 

In this work we shall consider only the GPDs corresponding to a Dirac structure F which conserves the quark 
chirality, that is F = and F = ^75. The associated matrix elements may be parametrized in the following way: 




FIG. 1: "Handbag" diagram. 
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where ujv is a nucleon spinor and H, H , E are the twist-2 chirality even GPDs. In the forward limit for which 
^ = and q^ — they reduce to the ordinary parton distributions, namely the longitudinal momentum q{x) and the 
helicity A.q(x) distributions are given by: 



q{x) ^ H{x,0,0), and Aq{x) = H{x,0,0) . 
The first few Mellin moments of these parton distributions are of particular interest 

x"-^^q{x) dx 



J x""-^ Aq{x) dx . 



(4) 
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Since as already mentioned, matrix elements of the light-cone operator as defined in Eq.([T]) cannot be extracted from 
correlators in euclidean lattice QCD, the usual method is to proceed with an operator product expansion of this 
operator that leads to a tower of local operators given by: 



(7) 

(8) 



The curly brackets represent a symmetrization over indices and subtraction of traces. The computation of the 
matrix elements of these operators on the Euclidean lattice can be done with standard techniques. The case n = 1 
amounts to calculating the elastic form factors of the vector and axial- vector currents and the results are reported 
in Ref. 0, In this work we concentrate on the n = 2 moments, i.e. matrix elements of operators with a single 
derivative. The matrix elements of the these operators are parametrized in terms of the generalized form factors 
(GFFs) ^20(9^), B2o{q^), C2o((?^) andi2o(g^), i32o(g^), according to 



{N{p',s')\0^l'\N{p,s)) = un{p',s') A2o{q')j^''P-^ +B2o{q') 
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Note that the GFFs depend only on the squared momentum transfer — {p' ~p)^ which imphes that the moments 
of the GPDs are polynomial in ^. In the forward limit we have v42o(0) — {x)q and ^20(0) = {x)Aq, which are 
respectively the first moment of the unpolarized and polarized quark distributions. Knowing the GFFs one can 
evaluate the quark contribution to the nucleon spin using Ji's sum rule: J« = i[A^o(0) + B|o(0)]. Moreover, usmg 
the measured or calculated value of the quark helicity AE^ — g\ the decomposition = + allows to study 

the role of the quark orbital angular momentum L'^ . 



II. LATTICE EVALUATION 

Twisted mass fermions [l^ provide an attractive formulation of lattice QCD that allows for automatic 0{a) im- 
provement, infrared regularization of small eigenvalues and fast dynamical simulations (l9j . For the calculation of the 
moments of GPDs, which is the main focus of this work, the automatic 0{a) improvement is particularly relevant 
since it is achieved by tuning only one parameter in the action, requiring no further improvements on the operator 
level. 

The action for two degenerate flavors of quarks in twisted mass QCD is given by 

^ = ^9 + ^ [fiv+mcrit+nsT^M] X{x) , (10) 

X 

where is the Wilson Dirac operator and we use the tree-level Symanzik improved gauge action Sg [20j|. The 
quark fields x in the so-called "twisted basis" obtained from the "physical basis" at maximal twist by a simple 
transformation: 

V'=-^[1 + »A5]X and V;=x-^[l + iT^75]- (11) 

We note that, in the continuum, this action is equivalent to the standard QCD action. As we pointed out, a crucial 
advantage is the fact that by tuning a single parameter, namely the bare untwisted quark mass to its critical value 
TOcriti a wide class of physical observables are automatically 0{a) improved [IE, E^, IMI . A disadvantage is the explicit 
flavor symmetry breaking. In a recent paper we have checked that this breaking is small for the baryon observables 
under consideration in this work and for the lattice spacings that we use (23 - [26j . To extract the GFFs without needing 
to evaluate the disconnected contributions we evaluate the nucleon matrix elements corresponding to the operators 
defined by 



where from now on we use the notation Oya'^" and O^a'^" to denote the vector and axial- vector operators with flavor 
index a. These matrix elements receive contributions only from the connected diagram for a — \, 2 and up to C(a^) 
for a = 3 prj . In particular, we consider the isovector combination with a = 3 for which the form of the operators 
remain the same in the physical and twisted basis. In order to find the spin carried by each quark in the nucleon 
we also analyse the isoscalar one-derivative vector and axial-vector operators. The latter receive contributions from 
disconnected fermion loops, which we neglect in this analysis. Simulations including a dynamical strange quark are 
also available within the twisted mass formulation. Comparison of the nucleon mass obtained with two dynamical 
flavors and the nucleon mass including a dynamical strange quark has shown negligible dependence on the dynamical 
strange quark ^28i] . We therefore expect the results on the nucleon moments to show little sensitivity on a dynamical 
strange quark. This is also confirmed by comparing our results to those where a dynamical strange quark is included. 

In this work we consider simulations at three values of the coupling constant with lattice spacings 0.056 fm, 0.07 fm 
and to 0.089 fm determined from the nucleon mass. This enables us to obtain results in the continuum limit. We also 
examine finite size effects by comparing results on two lattices of spatial length L = 2.1 fm and L = 2.8 fm [29l - [3l| . 
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{Xf,tf) ^ — ■ {Xi,ti 



FIG. 2: Connected nucleon three-point function. 



Correlation functions 



The GFFs are extracted from diniensionless ratios of correlation functions. The two-point and three-point functions 
are defined by 



G{q,tf - U) =^ e-*(-^/--^')-^rO„ {Utf,Xf)JpiU,x.)) (13) 



x.x 



where we consider kinematics for which the final momentum = 0. We drop tf — ti from the argument of the 
three-point function since it will be kept fixed in our approach. The projection matrices and F*^ are given by: 

r° = i(l+7o), F'= = iF°757fc fc= 1,2,3. (15) 
The proton interpolating field written in the twisted basis at maximal twist is given by 



1 



J(x)-^[1+Z75]e'"'' u'^' {x)C-i^d\x) u^{x). (16) 

where C is the charge conjugation matrix. We use Gaussian smeared quark fields [s^, [s^ to increase the overlap with 
the proton state and decrease overlap with excited states. The smeared interpolating fields are given by: 

V 

F = {t + aHY, 

3 

We also apply APE-smearing to the gauge fields C/^ entering the hopping matrix H. The parameters for the Gaussian 
smearing a and n are optimized using the nucleon mass as described in Ref. [6j]. 

For correlators containing the isovector operators the disconnected diagrams are zero up to lattice artifacts, and can 
be safely neglected as we approach the continuum limit. The detailed investigation of volume and cut-off effects will 
be performed on isovector quantities, for which no contributions are neglected. They can be calculated by evaluating 
the connected diagram, shown schematically in Fig. [2] A standard approach to calculate the connected three-point 
function is using sequential inversions through the sink [s^. The creation operator is taken at at a fixed position 
Xi—0 (source). The annihilation operator at a later time tf (sink) carries momentum p'=0. The current couples to a 
quark at an intermediate time t and carries momentum q. Translation invariance enforces q = — p for our kinematics. 
At a fixed source-sink time separation we obtain results for all possible momentum transfers and insertion times as 
well as for any operator 0^'^^"''^"\ with one set of sequential inversions per choice of the sink. We perform separate 
inversions for each one of the four projection matrices F** given in Eq. (jl5p . 

Using the two- and three-point functions of Eqs. (jl3p -(ll5 p and considering only one derivative operators we form 
the ratio 



i?-(F\,-,.) . ^^^X^ Gip tf-t)GiO,t)GiO,tf) 
G{0,tf) \ G{0,tf-t)G{p,t)G{p,tf)' 
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which is optimized because it does not contain potentially noisy two-point functions at large separations and because 
correlations between its different factors reduce the statistical noise. For sufficiently large separations tf — t and t — ti 
this ratio becomes time-independent (plateau region): 



lim lim 



(19) 



From the plateau values of the renormalized asymptotic ratio Il(r'^,q)]i — ZIl{T'^,q) the generalized form factors 
can be extracted. The equations relating Il{T^,q) to the GFFs are given in Appendix A. All values of q resulting in 
the same q^, the four choices of F and the ten orientations /i, v of the operator lead to an over-constrained system of 
equations which is solved in the least-squares sense via a singular value decomposition of the coefficient matrix. All 
quantities will be given in Euclidean space with ^ —q^ the Euclidean momentum transfer squared. The coefficients 
follow from the matrix-element decomposition given in Eq. Q and may depend on the energy and mass of the nucleon 
as well as on the initial spatial momentum p — ~q (see Appendix A). It turns out that both the operators with \x — v 
and n ^ y are necessary to obtain all three one-derivative vector form factors. Since those two classes of operators on 
a lattice renormalize differently from each other [3^1 , renormalization has to be carried out already on the level of the 
ratios. Although the one-derivative axial form factors can be extracted using only correlators with fi ^ ly, we use all 
combinations of /i, in order to decrease the statistical error. In Fig. Owe show representative plateaus for different 
momentum and indices fj, and v. 
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FIG. 3: Ratios for the one derivative vector (left) and axial vector (right) operator for a few exemplary choices of the momentum. 
The solid lines with the bands indicate the fitted plateau values with their jackknife errors. From top to bottom the momentum 
takes values p= (0,0,0); (0, 1,0), (1,0,0) and (1, 1,0). 

Since we use sequential inversions through the sink we need to fix the sink-source separation. The statistical errors 
on the three-point function are kept as small as possible by using the smallest value for the sink-source time separation 
that still ensures that the excited state contributions are sufficiently suppressed. We have tested different values of the 
sink-source time separation and for the final results we use the following values which correspond to t f — ti ^ 1 
fm 

/? = 3.9 : (tf - U)/a=l2 , ^ = 4.05 : (tf - ti)/a=l6 , ^ = 4.20 : (t/ - U)/a=l8. 



B. Simulation details 



The input parameters of the calculation, namely P, L/a and a/i are summarized in Table HI The lattice spacing a 
is determined from the nucleon mass, and the reader is referred to Ref. 0, [l3| for more details. Here we present only 
the final values, which are 

a;3=3.9 = 0.089(1)(5), a0=4.o5 = 0.070(1)(4), a^=4.2 = 0.056(2)(3) , 

where the first error is statistical and the second systematic. The pion mass values, spanning a range from 260 MeV 
to 470 MeV, are taken from Ref. [s^. At ~ 300 MeV and /3=3.9 we have simulations for lattices of spatial size 
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L=2.1 fm and L—2.8 fm allowing to investigate finite size effects. Finite lattice spacing effects are studied using three 
sets of results at (3—3.9, (3—A.05 and /3=4.2 for the lowest and largest pion mass available in this work. These sets of 
gauge ensembles allow us to estimate all the systematic errors in order to produce reliable predictions for the nucleon 
one-derivative GFFs. 



13 = 3.9, a = 0.089(1)(5) fm, r„/a = 5.22(2) 


24^ X 48, L = 2.1 fm 
32^ X 64, L = 2.8 fm 


afj, 0.0040 0.0064 0.0085 0.010 
No. of confs 943 553 365 477 
m„ (GeV) 0.3032(16) 0.3770(9) 0.4319(12) 0.4675(12) 
m^L 3.27 4.06 4.66 5.04 
afi 0.003 0.004 
No. of confs 667 351 

(GeV) 0.2600(9) 0.2978(6) 
m^L 3.74 4.28 


13 = 4.05, a = 0.070(1)(4) fm, ro/a = 6.61(3) 


32^ X 64, L = 2.13 fm 


a/i 0.0030 0.0060 0.0080 
No. of confs 447 326 419 
m„ (GeV) 0.2925(18) 0.4035(18) 0.4653(15) 
m^L 3.32 4.58 5.28 


/3 = 4.2, a = 0.056(1)(4) fm ro/a = 8.31 


32=* X 64, L = 2.39 fm 
48^ X 96, L = 2.39 fm 


a^i 0.0065 
No. of confs 357 

(GeV) 0.4698(18) 
rriTrL 4.24 
afi 0.002 
No. of confs 245 

(GeV) 0.2622(11) 
rrinL 3.55 



TABLE I: Input parameters {j3,L,afi) of our lattice calculation and corresponding lattice spacing (a) and pion mass (m.^). 



C. Renormalization 

We determine the renormalization constants for the one-derivative operators non-perturbatively, in the RF-MOM 
scheme [sT]. We employ a momentum source [ssj and perform a perturbative subtraction of 0{a^) terms [stI. [39|. 
This subtracts the leading cut-off effects yielding only a very weak dependence of the renormalization factors on (ap)^ 
for which the (ap)^ — > limit can be reliably taken. It was also shown with high accuracy that the quark mass 
dependence is negligible for the aforementioned operators. We find the values 

Z!^='' = 0.970(26), 1.013(14), 1.097(6) 

Z^f" ^ 1.061(29), 1.131(18), 1.122(10) (20) 
Zif" = 1.076(1), 1.136(0), 1.165(10) 

at /3=3.9, 4.05 and 4.2 respectively. These are the values that we use in this work to renormalize the lattice matrix 
element. 

III. LATTICE RESULTS 

In this section we present our results on the nucleon generalized form factors ^2o(Q^), ^2o(Q^), C'2o(Q^) and 
^2o(Q^), B2o{Q^)- We examine their dependence on the lattice volume and spacing, as well as, on the pion mass. We 
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also compare with recent results from other collaborations. In particular, we discuss lattice artifacts for the results 
on the isovector combination for the renormalized nucleon matrix element of the one-derivative operators 

in the MS scheme at a scale n = 2 GeV. 

In order to obtain some estimates on the spin content of the nucleon we also analyse the isoscalar parts of the 
spin-independent and helicity quark distributions, which however neglect the disconnected contributions. 

The GFFs ^20 (Q^ = 0) and A2o{Q^ = 0) are computed directly from the matrix elements, whereas B2o{Q^ — 0), 
C2o(Q^ = 0) and B2o{Q^ — 0) are obtained by linearly extrapolating the ^ data. C20 is consistent to zero 
within error bars for all momentum transfers. 



A. Finite volume effects 



In order to access volume effects we compare in Fig. |4] results on the moments {x)u-d and (x)Ati-Ad computed 
on different lattice sizes as a function of m^. As already mentioned, both these quantities are directly obtained at 
= and require no assumption on their Q^-dependence. Alongside our results we also show results using Np = 2 
clover fermions [i^ (preliminary) , Np = 2-1-1 domain wall fermions (DWF) and domain wall valence quarks on 
an Np = 2 + 1 staggered sea (hybrid) [i^ . 
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FIG. 4: {x)u-d and |{a;)Ati-Ad using twisted mass fermions (this work), Np ~ 2 clover fermions (preliminary), hybrid [42 
and DWF 



The results shown in Fig. [S] using twisted mass fermions correspond to a pion mass of about 300 MeV and are 
computed on lattices of spatial L with LttItt = 3.3 and LttItt = 4.3. As can be seen, results on these two lattices for 
both {x)u-d and (a::)AM-Ad are consistent. The LHPC using a hybrid approach and ^ 350 MeV has very accurate 
results at two lattices with LttItt = 4.5 and Lm^ = 6.2. No volume effects are seen for both vector and axial-vector 
first moments. The QCDSF collaboration has preliminary results for {x)u-d using clover fermions for ~ 270 MeV 
with Lm-jr = 3.4 and Lm^r = 4.2, which are consistent. Finally the RBC-UKQCD results with domain wall fermions 
with Lm^r = 3.9 and Lm^ = 5.7 show no volume effects for both {x)u-d and {x)Au-Ad [lH. The conclusion that 
we draw from this comparison is that finite volume effects on {x)u-d and (cc)A«-Aci are insignificant at our current 
statistical precision for lattices that satisfy im^r > 3.3. 

In Figs. Owe compare results on the GFFs A2q{Q'^), B2o{Q^), A2a{Q'^) and B2o{Q^) using twisted mass fermions 
for ttItt ~ 300 MeV for our two spatial lattice sizes oi L = 2.1 fm and L ~ 2.8 fm. The lines shown are linear fits to 
the Q^-dependence. As can be seen, for both v42o(Q^) and A2q{Q^) one can not ascertain any volume dependence. 
For B2o{Q^) and i?2o('9^) the statistical errors are larger and the linear fits show larger spread with the change in the 
spatial volume. However, given the large statistical uncertainties, it is difficult to quantify any volume dependence. 
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FIG. 5: Left panel: The isovector GFFs A20 and ,820 for pion mass ~ 300 MeV for a lattice of spatial length L 
L — 2.8 fm. Right panel: The axial-vector GFFs A20 and B20 for L — 2.1 fm and L = 2.8 fm. 
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B. Cut-off effects 



In order to examine cut-off effects we compare in Fig. [5] our results obtained at tfie lowest and highest pion mass 
that we have considered in this work, namely ^ 260 MeV and m^r ~ 470 MeV. We show results on the quantities 
^2o(Q^) and A2o{Q^) since these have smaller statistical errors than B2o{Q^) and B{Q^). For the heavier mass, where 
we have results at all three lattice spacings, there is no visible dependence on the lattice spacing especially at low 
Q^-values. For the lightest pion mass of m^r = 260 MeV we have results at the largest and smallest lattice spacings. 
The results are in good agreement although some deviations are seen at larger Q^-values. Thus, within our current 
statistical errors, one may conclude that no significant cut-off effects are observed. 



C. Quark mass dependence 

The mass dependence of A20 and v42o is shown in Fig. [71 Although the dependence on the mass is weak in the range 
of pion masses spanned, the tendency is for the values of the GFFs to decrease with decreasing pion mass. There is 
also a tendency for an increase in the slope for both ^20 and as the pion mass decreases. 



D. Comparison with other discretization schemes 



In order to compare lattice data using different discretization schemes one would have to first extrapolate to 
the continuum limit. However, given that the cut-off effects are small for lattice spacings of about 0.1 fm, lattice 
results for different values of a using a number of improved discretizations can be directly compared. In Fig. [5] we 
show results on the spin-independent and helicity moments using twisted mass fcrmions, in the hybrid approach 
obtained by LHPC [42], with Np = 2 clover fermions by the QCDSF collaboration J^] (preliminary results) and with 
DWF by the RBC-UKQCD collaborations [4l[. There is good agreement among lattice results. The very accurate 
results obtained using a hybrid action of domain wall valence and Np = 2 + 1 staggered fermions [i^l tend to be 
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FIG. 6: Left panel: The GFFs A20 and A20 for our three lattice spacings at ~ 470 MeV. Right panel: A20 and A20 for 
a = 0.089 fm and a = 0.056 fm at ~ 260 MeV. 



lower compared to the other data. One difference between them and the other results presented is that they are 
perturbatively renormalized. It was shown in Ref. 4l| that perturbative renormalization can lead to lower values. 
The spread in the values of the lattice results is shown to be reduced by taking a renormalization free ratio leading 
to a better agreement among lattice data with Lm^ > 4 (43j . In particular, constructing a renormalization free ratio 
brought the hybrid data in agreement with our results using twisted mass fermions and those using clover fermions by 
QCDSF. Lattice values for {x)u-d — ^20 (Q'^ — 0) although c omp atible are higher from the phenomenological value 
(xu-d) ^ 0.16. The very recent preliminary result by QCDSF j4(l| at to^ ~ 170 MeV remains higher than experiment 
and highlights the need to understand such deviations. A similar conclusion holds for the helicity moment. 

In Figs. [3] we compare our results for the GFFs with pion mass 300 MeV with those obtained using a hybrid action 
by LHPC and clover fermions by QCDSF and pion mass of 355 MeV and 350 MeV, respectively. The results show 
an overall agreement, with the data of LHPC somewhat lower than the other two sets. Once more the fact that 
both our results and those of QCDSF are renormalized non-perturbatively, while those of LHPC are renormalized 
perturbatively, might explain this difference. Moreover, in our determination of the renormalization constants, we 
have subtracted ©(a^) terms perturbatively to reduce lattice artifacts [STj . 



IV. CHIRAL PERTURBATION THEORY 



In order to make a direct comparison with experiment we need to extrapolate to the physical point. We only 
perform chiral extrapolation of GFFs at = 0. We first perform this extrapolation using our lattice results directly, 
since, as we have discussed in the previous section, cut-off effects are small. In the next section we will perform a 
continuum extrapolation and verify that indeed the values we find at the physical point are compatible. 

Within heavy baryon chiral perturbation theory (HB^PT) [49l [soj the expressions for the mTr-dependence of A20 
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from Ref. ^] for the unpolarized and from Ref. [45l.l46| for the polarized first moment. 



and A20 are given by: 
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where we take = 1 GeV^. The best fit is shown in Fig. [TUl where the width of the band is computed through 
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= 350 MeV. 



a super-jackknife analysis [42] • As can be seen, the fits yield a value higher than experiment for both observables. 
LHPC carried out a combined chiral fit using OijP') covariant baryon chiral perturbation theory (CBxPT) [ET'l to 
^20 J -B20 and C20- The mass of the nucleon at the chiral limit is used as input to the fits. They obtained a value 
for A20 in agreement with experiment [43 |. In order to compare with their analysis we also perform a combined fit 
to ^20(0), B2o(0) and C2o(0) within CBxPT [U]. The CBxPT fits are shown by the bands in Fig. [m As can be 
seen, they also provide a good description to the lattice data but, in the case of j42o, CB^PT leads to an even higher 
value at the physical point. Therefore, the discrepancy between our lattice results and the experimental value is not 
resolved. In Appendix B we collect the formulae used for the chiral extrapolations. The actual renormalized lattice 
data are tabulated in the tables of Appendix C for the isovector GFFs. 



V. CONTINUUM EXTRAPOLATION 



In order to study the continuum extrapolation we use the simulations at our three lattice spacings at the smallest 
and largest pion mass. We first interpolate the GFFs at the three values of /3 to a given value of the pseudoscalar 
mass in units of rg. We take as reference pion masses the ones computed on the finest lattice and interpolate results 
at the other two /?- values to these two reference masses. 

As already mentioned, < x >u-d= ^20 (0) and {x)Au~Ad — ^20 (0) are calculated directly at Q-^ — requiring no 
fits. We therefore choose these quantities to examine their dependence on the lattice spacing since this choice avoids 
any systematic errors due to the extrapolation to = 0, which would require the adaptation of an Ansatz for the 
Q^-dependence. 

Having determined the values at a given reference pion mass, we perform a fit to these data using the form 
y{a) = y(0) + c(a/ro)^. The resulting fits are shown in Fig. [T^l Setting c = we obtained the constant line also 
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The 



shown in the figure. As can be seen, for both large and small pion masses allowing a non-zero slope yields a value in 
the continuum limit that is in agreement with that obtained using a constant fit. This analysis shows that finite a 
effects are small for both large and small pion masses and extrapolation to the continuum limit using a constant fit is 
acceptable. For the intermediate pion masses we therefore obtained the values in the continuum by fitting our data 
at /3 = 3.9 and /? = 4.05 to a constant. For comparison we also perform a similar analysis for B2o{0) which requires 
fitting the Q^-dependence. The qualitative behavior is similar to the one observed for ^20 (0) and i?2o(0). 

Having results at the continuum limit we perform a chiral fit using HBxPT. The resulting curves are shown in 
Fig. [inland still produce a value at the physical point that is higher than the experimental value. In fact, the value 
obtained at the physical point for both vector and axial-vector moments is in agreement to the one extracted using 
the raw lattice data. This provides an a posteriori justification of using continuum chiral perturbation theory directly 
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FIG. 13: {x)u-d (left) and {x}Au-Ad (right) extrapolated to a = as a function of m,r. The blue band is the chiral fit using 
HBxPT.The physical point, shown by the asterisk, is from Ref. [U for the unpolarized and from Refs. [45l.[4^ for the polarized 
first moment. 



on the lattice data obtained at our three lattice spacings to perform the extrapolation to the physical point in the 
previous section. 



VI. PROTON SPIN 



In order to extract information on the spin content of the nucleon one needs to evaluate the isoscalar moments 
^20^'* and B2q'^ since the total spin of a quark in the nucleon is given by 

J'^liAUO) + Blom- (23) 
The total spin can be further decomposed into its orbital angular momentum and its spin component AS"^ as 

J« = IaS^ + L« (24) 

The spin carried by the u- and d- quarks is determined using AI]"+'' — A^q"^. In order to evaluate the isoscalar 
quantities one would need the disconnected contributions. These are notoriously difficult to calculate and they are 
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neglected in most current evaluations of GFFs. Under the assumption that these are small we may extract the 
information on the nucleon spin. 
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FIG. 14: Chiral extrapolation using CBxPT (left) and HB^PT (right) for the total spin carried by the u-and d- quarks. 




FIG. 15: Chiral extrapolation using HBxPT for the angular momentum and spin carried by the u-and d- quarks. The physical 
points, shown by the asterisks are from the HERMES 2007 analysis ,52]. 

In Fig.[Tl]we show our results for the isoscalar ^20(0)"^'': .620(0)"+'^ and C2o(0)"~'''^. Since as shown in the previous 
section, cut-off effects are small, we here perform a chiral extrapolation directly on the lattice data. Having both 
isoscalar and isovector quantities we can extract the spin J"" and J'^ carried by the u- and d- quarks. The results are 
shown in Fig.[T4l We show the extrapolation using both HBxPT and CBxPT both of which have the same qualitative 
behavior. As can be seen, the contribution to the spin from the d-quark is much smaller than that of the u-quarks. 
These results are in qualitative agreement with the recent results obtained using a hybrid action [i^. In Fig. (TS] 
we show separately the orbital angular momentum and spin carried by the u- and d- quarks. Both our results and 
those of LHPC [421 qualitative agreement as far as the spin is concerned. For the orbital angular momentum 

we obtain higher values for both d- and u-quark (less negative). Thus we obtain a total positive compared to 

a small negative value in the case of LHPC. After chiral extrapolation, the value obtained at the physical point is 
consistent with zero in agreement with the value determined by LHPC. For the spin contribution AI]"+'' our value 
is lower at the physical point as compared to that obtained by LHPC. We summarize the values for the total spin, 
orbital angular momentum and spin in the proton at the physical point in Table HIl 



15 





CBxPT 


HBxPT 


PYT^PT"! TTI PTii" 
CAJJCJ. illiCliLj 


ru — d 
U 








ju + d 
u 


















-fi ni ^("s"! 






AE""'' /2 












0.148(5) 


0.208(9) 


AE"/2 




0.305(7) 


0.421(6) 


AE'*/2 




-0.157(5) 


-0.214(6) 






-0.258(5) 




j^U'\-d 




-0.025(53) 








-0.141(30) 








0.116(27) 





TABLE II: Values of nucleon spin observables at the physical point using CB^PT and HB^PT and from experiment [53 |. 



VII. CONCLUSIONS 

We have performed an analysis on the generalized form factors, ^2o(Q^), S2o(<3^), C2o{Q^), ^2o(Q^), ^2o(Q^), 
extracted from the nucleon matrix elements of the one-derivative vector and axial-vector operators using two degen- 
erate flavors of twisted mass fermions. Our results are non-perturbatively renormalized and they are presented in the 
MS scheme at a scale of 2 GeV. To investigate volume and cut-off effects we have used the isovector combinations, 
which can be calculated without the necessity to evaluate disconnected contributions. Our main conclusion regarding 
cut-off effects is that they are small within the current accuracy of about 5-10% and for lattice spacings smaller than 
0.1 fm. Similarly, no systematic volume effects are seen. Given the small cut-off effects one can compare lattice results 
directly using different discretization schemes. The comparison of the results using Np = 2 twisted mass fermions 
with the results obtained using Np = 2 clover fermions by the QCDSF [40|] shows agreement. Both the results of this 
work as well as those by QCDSF are non-perturbatively renormalized. We also compared our results with Np = 2 + 1 
domain wall fermions ^l"!. Again there is agreement without any indication of any systematic effect from including 
a dynamical strange quark. Our results at three values of the lattice spacing allow for a continuum extrapolation. 
By interpolating results to a reference mass in units of ro and performing a linear extrapolation in it was shown 
that the values obtained are consistent with those obtained with a constant extrapolation. This has been verified for 
both the heaviest and lightest masses used in this work. Furthermore, if one performs chiral fits to the extrapolated 
continuum results one finds a value at the physical point consistent with the one obtained using directly the lattice 
data at finite lattice spacing. This provides a consistency check that cut-off effects for a lattice spacing less than 
0.1 fm are smaller than our current statistical errors. 

Having established that both volume and cut-off effects are small for the isovector quantities for which only con- 
nected contributions are needed, we analyse the corresponding isoscalar quantities using directly our lattice data. Of 
particular interest here is the spin content of the nucleon. The disconnected contributions to the isoscalar quantities 
are not included. We find that the spin carried by the d-quark is almost zero whereas the u-quarks carry about 50% 
of the nucleon's spin. This result is consistent with other lattice calculations p^ . 

For the chiral extrapolations of these quantities we use HB^PT and CB^PT theory. In both cases, our results on 
the momentum fraction and hclicity moment at the physical point are higher than their experimental value. Such 
discrepancies are also observed in the case of the nucleon axial charge and they need to be further investigated. 
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Appendix A: Expressions for the extraction of GFFs from lattice measurements 

Wc collect here the expressions relating the plateau values to the GFFs A20, -B20, C20 and ^20, -820- The index V 
(A) refers to the vector (axial-vector) one-derivative operator. All relations are given in Euclidean space. 
Vector: 



3En El mN\ , r, ( En E% El itinX 
- 7^ - ^ + ^20 O - — + —2- + — — - - — 



\ 8 4mjv 8 y V o 8m^ mmN 16 / 

, ^ r ( En El El TUN \ , s 

no?(r",g-) = 0, (26) 

\ 8 8 4:mN J V 16mAr 16 8ml 8mjv/ 

+ C20C + ^ + + , (27) 



4mjv 4 2m|f 2 tun J ' 

n^'=(r", g) = A20 c (z IhJipi^) +B20C (t '^-^pi^ip ') , (28) 

V 4mAr / V 4mjv 



Ily{T ,q) = A20C [-1 — - I — — B20 C -« + « -^—2- + C20 C [t^-t ir-T- . (29 
V 4 4 mjv / V 8 8 / V 2 2 ml ) 

= M.C (.„„ (I + I^)) +i>»c (I + I^)) . (30) 

n«,r.„T = .,„c|i| + B.c(-«^ + ii5)+c»c(^ + |i5) , (31) 

n«(r",J| = A„ C fi < "' "l + B^o C (i i V (32) 

\ 8mjv 8mjv / V 8mjv 8mjv / 

Axial - Vector : (33) 

n^^(ro,0 = 0, (34) 

n-(r», = i2o c f ^„ + -h ^ V * 1^1 + ^ f ^ . (35) 



4 ' 8mjv ' 8 ) ' 8mN J ' V 



nl^(r",5-) = i.oC f<5„, ff + (f + 1^)) +^2oC (-^) , (36) 

V \8 8tojv / V8 8mN J J \ 8m^ / 

where C = -^2ml/ {En{Em + m^)), El = ml + (f and the Latin indices k,j,n denote spatial directions 1, 2, 3 and 
k ^ j. A summation is implied over the index p. 
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Appendix B: Chiral perturbation theory results 

For convenience we collect in tliis appendix the results of HB^PT results taken from Ref. [s^] for the isovector 
(/ — 1) and isoscalar (/ ~ 0) first moments and axial charge: 
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We note that the expressions for A2q^ and A2q^ are the same as those given in Eq. (up to a redefinition of C, 
(7 and cg, cg). We have included them here using the notation of Ref. [ss'] for completeness. 
We performed a combined fit to the following CBxPT results taken from Ref. [5l|: 
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19 



^20 (0) - a2o + 4m,^2 ^^^^^^^ 
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where Mq is the mass of the nucleon at the chiral hmit. 
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Appendix C: Numerical results for the isovector sector 
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TABLE III: Results on A20 and B20 form factors at /3 = 3.9, lattice size: 24^ X 48 
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TABLE IV: Results on A20 and B20 form factors at /3 = 3.9, 



lattice size: 
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TABLE V: Results on A20 and B20 form 



factors at 
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TABLE VI: Results on A20 and B20 form factors at /3 = 4.2 



